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Using the Lagrangian formalism, we solve analytically the equations of motion for current-induced 
domain-wall dynamics in a ferromagnet with Rashba spin-orbit coupling. An exact solution for the 
domain wall velocity is provided, including the effect of non-equilibrium conduction electron spin- 
density, Gilbert damping, and the Rashba interaction parameter. We demonstrate explicitly that 
the influence of spin-orbit interaction can be qualitatively different from the role of non-adiabatic 
spin-torque in the sense that the former is sensitive to the chirality of the domain wall whereas 
the latter is not: the domain wall velocity shows a reentrant behavior upon changing the chirality 
of the domain wall. This could be used to experimentally distinguish between the spin-orbit and 
non-adiabatic contribution to the wall speed. A quantitative estimate for the attainable domain 
wall velocity is given, based on an experimentally relevant set of parameters for the system. 

PACS numbers: 



I. INTRODUCTION 

The study of domain- wall motion in ferromagnetic ma- 
terials has attracted much interest in recent years. Be- 
sides its allure from a fundamental physics viewpoint, 
electric control of magnetic textures is attractive in terms 
of potential new applications such as magnetic memory. 
A key concept in the context of domain wall motion is 
the so-called spin-transfer torque [P-^ : in essence, it con- 
sists of a transfer of a transverse spin-current component 
to the ferromagnetic order parameter which may occur 
in a non-coUinear magnetization configuration. Active 
control over domain-wall motion is a chief objective in 
terms of realizing the "magnetic race-track" technology 
put forward in . Other ways to manipulate domain wall 
motion include energy redistribution in the presence of 
an external magnetic field, by applying microwave radia- 
tion [^^, and by using magnons ^10H13] . Moreover, the 
concept of spin-transfer torque has recently been studied 
in antiferromagnets p!4lfl9] including the possibility to 
move domain- walls. 

Interestingly, it turns out that spin-orbit interactions 
can substantially modify both the spin-transfer torque 
and the resulting domain wall motion [201 HIl [IS] . The 
combined influence of a magnetic exchange field together 
with spin-orbit interaction, typically taken in the Rashba 
form, can be shown to give rise to a non-equilibrium 
spin density perpendicular to the injected current flow. 
In turn, this gives rise to an effective magnetic field 
which causes magnetization dynamics and, for suffi- 
ciently strong current density, magnetization switching. 
The effective " spin-orbit torque" arising in this manner is 
qualitatively different from the conventional spin-transfer 
torque due to the different mechanism at hand: it does 
not require the presence of non-coUinear magnetic ele- 
ments and will cause magnetization dynamics in a single 
ferromagnetic layer [28]. In addition, it is important to 
note that whereas different types of domain walls behave 
in the same way in the absence of spin-orbit interactions, 
the exact magnetization texture plays a key role when 



spin-orbit coupling is present due to the coupling between 
the electron motion and the spin torque. 

In light of the above, several experimental and theoreti- 
cal works has recently explored the influence of spin-orbit 
interactions on magnetization dynamics in various mag- 
netic systems p2l430j . Since the addition of spin-orbit 
terms in the equations of motion for the magnetization 
texture complicates their solution, the large majority of 
these works have relied on numerical methods to solve 
the Landau-Lifshitz-Gilbert (LLG) equation. In this pa- 
per, we utilize the Lagrangian formalism in order to write 
down and solve analytically the equations of motion for 
a domain wall within a collective-coordinate description. 
The analytical nature of this approach allows us to iden- 
tify a transparent expression for the domain wall velocity 
and how it depends on parameters such as the spin-orbit 
interaction, exchange field, and Gilbert damping of the 
system. Alternatively, one could have derived this result 
via the LLG equation, but the present formalism makes 
it easier to accommodate non-equilibrium spin-density 
terms and higher order corrections due to the spin-orbit 
interaction. We show from the analytical expression that 
the presence of spin-orbit interactions renders the domain 
wall velocity to be chirality sensitive, in effect depend- 
ing on whether the wall changes magnetization direction 
from positive to negative or vice versa along the direction 
of the current. We provide an estimate for the magnitude 
of the domain wall velocity using a set of experimentally 
relevant parameters and show that the velocity behaves 
qualitatively differently depending on the chirality of the 
domain wall and displays reentrant behavior. Our finding 
suggests a way to experimentally distinguish between the 
non-adiabatic and spin-orbit contribution to the domain- 
wall speed. 



II. THEORY 

We consider first a ferromagnetic domain wall with an 
easy (hard) axis of magnetic anisotropy |36j along the 
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FIG. 1: (Color online) Sketch of the magnetization textures 
considered in this paper. The domain wall nanowire extends 
along the a;-axis. (a) Hard axis along current direction (x) and 
easy axis along z-direction. In equilibrium, the magnetization 
of the domain wall rotates in the j/z-plane exactly as in Ref. 
|21| and similarly to Ref. Out of equilibrium, i.e. under 

a current bias, a finite component may be acquired along the 
a;-axis as indicated by the tilt angle (j>- This geometry is the 
main focus of this article, relevant for nanowires with perpen- 
dicular magnetic anisotropy. (b) Hard axis along j/-direction 
and easy axis along z-direction. In equilibrium, the magne- 
tization rotates in the a;2:-plane. Out of equilibrium, a finite 
component may be acquired along the y-axis as indicated by 
the tilt angle 0. 

y-axis (x-axis). An electric current is injected along the 
X-axis [see Fig. [l^)], and the inversion symmetry is as- 
sumed broken in the z-direction. This gives rise to a 
Rashba spin-orbit coupling term, which influences the 
magnetization dynamics. More precisely, it can be shown 
[20] that the spin-orbit coupling generates an effective 
magnetic field 

-H'SOC OC aRZ X je (1) 

where is the Rashba interaction strength and j'e is 
the current density vector. In order to treat the domain 
wall as rigid in a collective-coordinate framework, thus 
making other modes of deformation apart from the tilt 
angle (j) irrelevant, the easy axis anisotropy energy K is 
assumed larger than its hard axis equivalent K± |31| . i.e. 
\K\ ^ \K±\. As is normally done, we do not take into 
account the effect of the end-boundaries of the nanowire, 
assuming thus that the domain-wall center is located suf- 
ficiently far away from these during its propagation. 
The starting point is the Lagrangian for such a Bloch 



domain- wall which was derived in Ref. . For this type 
of domain wall, the magnetization rotates in the plane 
perpendicular to the extension of the wire (and current 
direction) , similarly to Ref. [55] (the only difference from 
Ref. [21] is which of the perpendicular axes that is the 
easy one, in both cases the hard axis is along the current 
direction). It reads: 

£ = (# - sin^ 0) - 2XxJ- - s^'xt/- - 0jf — + TiX)) 

jJi \ fJL J 

(2) 

with the definition 

Above, 4) is the tilt angle of the domain wall (see Fig. 
[T]) while X is the normalized position of the center of 
the domain wall {x = X/L where L is the wall thick- 
ness). The quantities A and J are the Rashba interac- 
tion strength and the current density normalized against 
niL/h? and evc, respectively, with Vc being the drift ve- 
locity of electrons at the intrinsic threshold current for 
A//i = 1 without Rashba interactions. Moreover, A is 
the exchange splitting, /i is the Fermi energy while m 
is the electron mass. Finally, represents the constant 
non-equilibrium spin density induced by the applied elec- 
tric field generating the current. 

III. RESULTS AND DISCUSSION 

We will model dissipation in this system by a Rayleigh 
dissipation function of the form [3T] 

W=f(i' + '^')- (4) 
The Lagrange equations are obtained via 

d dc dc dw . 

d^a^-a^^-^''^^^"''^^- 

where t is dimensionless time-coordinate normalized 
against l/vc- Defining for convenience 

c=-+T{X) (6) 
M 

we obtain the following Lagrange-equations, which were 
studied numerically in [21] : 

A 

(j) + ax — s^(j> — 2XJ — , X — acf) = —s^x + sin 2<j> + cJ. 

(7) 

The role of the dissipative spin-transfer torque (also 
known as non-adiabatic or /3-torque) will be addressed 
later on - we will see that it plays a similar role as the 
spin-orbit coupling, but with one important difference. 
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FIG. 2: (Color online) Plot of the terminal domain wall ve- 
locity using a = 0.005, L = 75 nm, m = 0.04me, \i = 0.05 
eV, A = 0.02 eV, and Vc = 150 m/s. (a) Positive chirality 
A > 0. (b) Negative chirality A < 0. 



By combining the above equations, we are able to elimi- 
nate the i-dependence and obtain a first-order differen- 
tial equation for the tilt angle: 



a 



sin 20 



cJ 

1 + sv 



(8) 



This equation may be cast into integral form as follows: 



a + sin 24> 



t. (9) 



We define the quantities 



a = cJ 



2AA J(l + sy) 

b^-a-a-^[l~{syf] 



(10) 



The formal solution of this integral is obtained after some 
algebraic manipulation: 



tan </) = -a"^ -I- a"^ - 1 tan(6 ^t\/a^ - 1 + Cq), 

(11) 

where Co is an integration constant to be determined 
from the initial conditions. In particular, (j>{t = 0) = 
and (j){t = 0) = tt/2 yield Cq = atan(l/-\/a2 — 1) and 
Co = 7r/2, respectively. 

Having obtained an explicit expression for the tilt 
angle, we are now in a position to identify the time- 
dependence of the domain- wall center x, thus also ob- 
taining the domain -wa ll velocity x. To do so, we first 
obtain (j> from Eq. ( 11 ) as: 



- 1) 



'■{b-HVa^+Co) 



(12) 



ah 1 + [-a-i +a-Va2 - ltan(fe-i<%/a2^ + Co)]2' 
Substituting this into the first Lagrangc-equation, one obtains an explicit expression for the domain-wall velocity vuw 



cos-^ {b-HVa^ - I + Co) 



AA ,, (a^~l) 

'^DW = (1 — S'') ; , , . 

a/i aba 1 + [-a-'^ + a-^VaF^Un{b-HVa?^ + CoW 



(13) 



Eqs. (11) and (13) determine an exact analytical expression for the time-dependence of the domain wall tilt-angle 



and velocity, respectively, which we will analyze in more det ail below. 



r 



By integration, Eq. ( 13 ) may be used to identify the taneous position of the domain-wall center 



time-evolution of the domain-wall center: 

„ A A ii-sy) f 

X = -2Jt ^ '- / dt^. 

a II CK J 



(14) 



The integral over (j) is evaluated by making use of the 
formula 



dycos ^(j/ + 7) 
[l + (77 + Ctan(a; + 7)) 



2 =C ^atan[Ctan(a; + 7) + 77], 

(15) 



resulting in the following equation describing the instan- 



X — xo — —2Jt- 



A A 



1 - s'^ 



- X 



an a 
tan-i[a-Va^ - 1 tan(&"4v/a2 - 1 + Cq) 



(16) 



with xq being an integration constant related to the ini- 
tial position of the domain-wall. 

It is worth noting here that there exists a particularly 
simple solution to the equation set Eq. ([t]) in the spe- 
cial case of a constant tilt angle, i.e. a domain-wall that 
preserves its shape and magnetization direction and thus 
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only propagates. This amounts to setting cf) = 
dictates that the tilt angle must be constant: 

sin(20o) = -cJ- 2(1 + 5'') J- — , 

a fi 

and gives for the domain- wall velocity: 



0, which 



(17) 



vdw 



-2J- 



a fi 



(18) 



This is identical to the first term in Eq. ( 13 1. Physically, 



this implies a constant drift velocity of the domain-wall 
under the application of a current. Interestingly, it is 
seen that wdw = in the absence of A, meaning that the 
spin-orbit interaction is fully responsible for the domain- 
wall motion. In a more general scenario where the tilt 
angle is not restricted to being constant, i.e. allowing for 
domain-wall deformation, the general expression for the 
velocity is given by Eq. (13 1. As we shall discuss later. 



this corresponds to the Walker breakdown threshold. 

The necessity of spin-orbit coupling to drive the 
domain-wall motion is a feature pertaining specifically 
to the Bloch-domain wall with anisotropy axis chosen 
as described in the beginning of this section [see Fig. 
[l^)]. In fact, for different choices of anisotropy directions, 
spin-orbit coupling mainly contributes as a quantitative 
correction to the domain-wall velocity without being a 
prerequisite for its existence. To see this, one may con- 
sider instead a domain-wall where the magnetic easy and 
hard axes lie along the z and y-axes, respectively [see 
Fig. [ija)]. In this case, the appropriate Lagrangian to 
consider is EH 



C = {(j)x — sin (p) + cos ( 
A 

— s'^xcj) — (j)J( — 



2Ami2 



{x cos (f) — sin 00) 
) + 7rAsin0JA/^. (19) 



Using the same approach as above, one arrives at the 
following set of integro-differential equations for the time- 
evolution of the domain-wall center and tilt angle: 



-ax 



(20) 



in addition to 



1 r 9 / T^S'^ 

ba ^[a +(1 + — sm 



— sin 20 — J( 



+ s^) V [ttA cos JA//Z 
)]-t-io. (21) 



^ 2AmL2 



Unlike the situation considered formerly, this equation 
set may be solved exactly analytically only under simpli- 
fying circumstances. For instance, by assuming that the 
final term in the denominator of the integral equation 
above dominates and moreover using s''' <i; 1, 7 = {y, 2}, 
one identifies the domain-wall velocity in the strong- 
current regime J ^ 1 as 



J 



'A 



(22) 



As seen, the presence of spin-orbit coupling in this case 
brings about a minor correction to the final velocity, es- 
pecially for a strong ferromagnet where A//i dominates 
the expression in parantheses. The domain-wall is driven 
directly by the current J with a velocity that increases 
with decreasing dissipation a. The above result for the 
domain wall velocity shows that the spin-orbit coupling 
has a qualitatively different effect on different domain 
wall textures. 



In the following, we focus on the more interesting case where the spin-orbit coupling influences qualitatively the 
domain- wall motion [Eqs. ( 11 1 and ( 13 )] and investigate the precise dynamics using the derived analytical expressions. 
For concreteness, we consider initia conditions such that the tilt angle of the domain-wall at t = is zero, meaning 



that we may write: 



tan = 



^'DW 



1 



\/ a? — 1 tan 
A A 



rWa? - 1 
V h 



at an 



1 



- 1 



-2J- 



a fj. 



a 



ab 



cos-2(6-ii\/a2 ~ 1 + atan(l/Va2 - 1)) 
1 + 4,[^/a^^tan(6-l<Va2 - 1 + atan(l/Va2 - 1)) - 1]2 ' 



(23) 



In the limit a — 00, we flnd that 0(t) — >■ and vuw ~^ 
as expected. The analytical expression for vuw above 
reveals that the velocity has non-monotonic behavior as 
a function of time. In particular, there is a resonance 
condition t — t^-cs at which the velocity increases in mag- 



nitude: 



(n+l/2)7r& b 
Va'^ - 1 - 



1' 



(24) 



assuming that a > 1. In effect, udw exhibits oscillations 
which persist even for its terminal behavior t — ^ 00. It is 
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therefore of interest to establish the average domain wall 
velocity by averaging over one period T = 7r6/Vt?^^: 

1 

(wDw) y y dtoDW- (25) 

Note that for a < 1, only the drift-term in udw survives 
as t — )• cxD [since cos^(±ii) — )• cxd in this limit]. Performing 



the integral with this in mind, one finds for arbitrary a 
that 

(fDw) = -2J— - ^— -^sign{a}7^e{v/a2 - 1}, (26) 
a/i ab 

which written out in terms of the original physical pa- 
rameters reads 



(^'Dw) = -2J 



AA (1 - s?')sign{ J} 



ajjL 



— — X sign<^ 



ajjL 



A 



X Tlci 



2AA(l + s2') 



A 



J-(A) 



J2 - 1 

(27) ' 



Eq. (27) is the main result of this work and constitutes 



a generally valid analytical expression for the terminal 
domain wall velocity taking into account both Rashba 
spin-orbit coupling, the non-equilibrium spin density, and 
Gilbert damping. For consistency, we have verified that 
a fully numerical solution of the equations of motion give 
identical results as the above analytical expression for the 
domain wall velocity. 

We will proceed to analyze this velocity quantitatively 
for a realistic set of parameters and investigate how it de- 
pends in particular on the applied current and the mag- 
nitude of the spin-orbit coupling. Before doing so, one 



should note that it follows from Eq. ( 27 ) that there exists 
both a threshold current Jr for which the second term in 



Eq. (27) is non-zero: 



\Jc\ 



2AA(1 



ajjL 



2\2 



mL^AJ 



(28) 



In the limiting case of zero spin-orbit interaction, one ob- 
tains \Jc\ — {A/fj,)^^ in agreement with previous studies. 
The presence of spin-orbit interaction is seen from the 
analytical expression of \Jc\ to reduce the threshold cur- 
rent monotonically with increasing A, consistently with 
the numerical study in Ref. |21j . This monotonic be- 
havior appears also when tuning the chemical potential 
fi: increasing fi lowers the polarization and increases the 
threshold current. 

Eq. ( 28 1 is in fact the Walker threshold value which 



separates the regimes of domain wall motion with a fixed 
profile, i.e. (j) = and the regime with a domain wall 
rotating its spatial profile as time increases, i.e. </> ^ 0. 
To see this, one may revert to the original equations of 
motion in Eq. ([7|. There exists a tilt angle which 
satisfies </> = if the following equation is satisfied: 



. /2XA{1 + sy) ^ 

sni2(j> = - ^ -J + cJ 



(29) 



Since the left-hand side varies between ±1, one may find 
a solution if the following equation holds: 



\J\ 



2AA(1 + gy) 



< 1. 



(30) 



which is completely equivalent to Eq. ( 28 1 after rewrit- 
ing. For larger currents J, there exists no time- 
independent solution (p and domain wall distortion (j) is 
now inevitable past the Walker breakdown. 

It has previously been suggested that the spin-orbit 
interaction and non-adiabatic spin-torque influence the 
magnetization dynamics in the same manner, since the 
latter may be included by substituting A — > /3 + A [21]. 
However, it was noted in Ref. that the chirality of the 
domain wall determines the effective sign of the spin-orbit 
coupling A in the equations of motion. Formally, this 
corresponds to taking the domain wall profile represented 
with polar angles (j> and 6 as M = Mq{— sm9 s'mipx + 
cos 9y + sin 9 cos (f)z) and performing the transformations 
(j) — >■ {—(p) and cos 6* {—cos 9). The parameter (p — 
(j){t) is the time-dependent tilt angle, whereas 9 is defined 
by 



sm9 — sech[(i — X(t))/LDw], 



(31) 



where x is the position along the magnetic wire, X{t) 
is the time-dependent center position of the domain- 
wall, whereas Ldw is the domain wall width. This sug- 
gests that the role of spin-orbit coupling is chirality- 
sensitive and in this regard differs qualitatively from non- 
adiabaticity. Below, we will investigate this effect and 
show that the chirality indeed gives rise to highly differ- 
ent behavior for the domain wall velocity. The chirality 
is changed in our analytical expression Eq. (27) by let- 
ting A -> (—A). We note that our conclusions remain 
unchanged even when including the /3-term, as long as it 
is small (which is typically the case). 

In the remainder of this paper, we fix the following pa- 
rameters: a = 0.005, L = 75 nm, m = O.OAme, fj, = 0.05 
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eV, A = 0.02 eV, and Vc = 150 m/s to model an exper- 
imentally realistic semiconductor system [32], where me 
is the bare electron mass. We restrict our attention to a 
scenario where the electron-spin density satisfies ^ 1. 
In Eq. (27), note that J and A are normalized quanti- 



ties. To make a quantitative estimate for the domain wall 
velocity, we define the non-normalized current — Jv^ 
and spin-orbit interaction A = hX/ (mL) which have units 
m/s and eV-m, respectively. Similarly, we restore the di- 
mension of the terminal domain wall velocity by defining 
(^Dw) — Vc{vdw) with units m/s. 

We show in Fig. [2] a plot of the terminal domain wall 
velocity as a function of the applied current for several 
values of the spin-orbit interactions. To illustrate the ef- 
fect of the chirality, we plot in (a) the case A > while 
in (b) A < 0. The latter results are consistent with the 
numerical study of Ref. [21], and shows that the spin- 
orbit interaction can greatly enhance the domain wall 
velocity at low currents up to the Walker breakdown. In 
the former case, however, the domain wall velocity be- 
haves differently when changing the current J^. The spin 
torque induced by the effective Rashba-field is now di- 
rected opposite to the conventional current-driven spin 
torque and there is a competition between the two. For 
non-zero A, the domain wall velocity starts by moving in 
the direction opposite to the current, whereas it eventu- 
ally changes sign with increasing J' (above the threshold 
current). Thus, experimentally observing a sign-reversal 
of the domain-wall velocity with applied current would 
be an indication of precisely this chirality sensitive spin- 
orbit coupling effect. We note that this sign-reversal of 
the wall velocity is different from the one predicted in 
[25], which originated from a Slonczewski-like spin-orbit 
torque proportional to /? and when considering a differ- 



ent type of domain wall profile. In our case, the field- 
like spin-orbit torque is sufficient to cause the velocity- 
reversal. 



IV. SUMMARY 

In conclusion, we have used the Lagrangian formalism 
to derive an exact analytical expression for the domain 
wall velocity in a spin-orbit coupled ferromagnet. We 
have shown that a chirality-sensitive domain wall 
velocity appears in this system, which qualitatively 
differs from the role of the non-adiabatic spin-transfer 
torque stemming from a spatial mistracking between the 
conduction electrons and the local magnetization. A 
candidate system for the observation of this effect would 
be hybrid structure comprised of a thin ferromagnetic 
wire in contact with a heavy metal and an oxide, which 
would break structural inversion symmetry and this pro- 
vide a gradient in the electric potential which generates 
the Rashba spin-orbit coupling. This type of structure 
has indeed recently been experimentally considered in 
Refs. [m US] (Pt/Co/AlOx) whereas similar systems 
were studied in Refs. [SSSS] (Ta/CoFeB/MgO). 
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